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ABSTRACT 

f(R) gravity, capable of driving the late-time acceleration of the universe, is emerging as 
a promising alternative to dark energy. Various f(R) gravity models have been intensively 
tested against probes of the expansion history, including type la supernovae (SNIa), the cos- 
mic microwave background (CMB) and baryon acoustic oscillations (BAO). In this paper we 
propose to use the statistical lens sample from Sloan Digital Sky Survey Quasar Lens Search 
Data Release 3 (SQLS DR3) to constrain f(R) gravity models. This sample can probe the 
expansion history up to z ~ 2.2, higher than what probed by current SNIa and BAO data. 
We adopt a typical parameterization of the form f(R) = R — aHn(-jk-)^ with a and (3 

constants. For /? = (ACDM), we obtain the best-fit value of the parameter a — —4.193, for 
which the 95% confidence interval that is [-4.633, -3.754]. This best-fit value of a corresponds 
to the matter density parameter fl m .o = 0.301, consistent with constraints from other probes. 
Allowing (3 to be free, the best-fit parameters are (a, (3) = (—3.777, 0.06195). Consequently, 
we give il m o — 0.285 and the deceleration parameter qo = —0.544. At the 95% confidence 
level, a and (3 are constrained to [-4.67, -2.89] and [-0.078, 0.202] respectively. Clearly, given 
the currently limited sample size, we can only constrain (3 within the accuracy of A/3 ~ 0.1 
and thus can not distinguish between ACDM and f(R) gravity with high significance, and 
actually, the former lies in the 68% confidence contour. We expect that the extension of the 
SQLS DR3 lens sample to the SDSS DR5 and SDSS-II will make constraints on the model 
more stringent. 
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1 INTRODUCTION 

In modern cosmology, one of the most striking discoveries is that 
our expanding universe is undergoing a phase of acceleration. 
The key observational results that support this disco very are; the 
luminosity-redshift relationship from SNIa su r veys dAstier et al.l 
l200d IPerlmutter et alj|"l999l: iRiess et al.lll998l [20041. l2007h. the 
CMB anisotropy spectrum dde Bernardis et alJl2000l : ISpergel et al.l 
120031 120071: Iffinshaw et ai] |2007h. the large scale structure from 
galaxy redshift surv eys dCole et all 120051; ISeliak et alj 120051 : 
iTegmark etail l2004ah and BAO feisenstein et al.l 120051) . With 
the combination of all these observational results, one gets the 
standard ACDM cosmology as: fix = 0, Qm = 0.27, = 
0.73. Amongst the matter content baryonic matter amounts to 
only 4%, other 23% is the so-called dark matter (DM), and the 
rest component of 73% dominate the universe that is often re- 
ferred to as dark energy (DE), which is a negative-pressure ideal 
fluid smoothly permeated the universe. It is worth noting that, 
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notwithstanding the standard cosmology is very well in agree- 
ment with the astrophysical data, the nature of the dark matter 
and dark energy remains mystery in the modern Cosmology and 
Physics, so that there are so many theoretical proposals to ac- 
count for them on the ground. Possible models t o explain the ac- 
celeration include a classical cosmo logy constant dWeinberg||l989 | : 
Carroll etalJI 1992b. Chaplygin gas jKamenshchik et al.ll200ll ; lzhu] 



2004l : IWu & Yul2007allbh , a wide variety of scalar-field models such 



as Quintess ence (e.g.. ICaldwell et alj] 19981: IWu & Zhandl2008h 



K-essence ( Armendariz-Picon et al 



2000, 



2001), Phantom field 



dCaldw ell 2002; Dabr owski et alj|2003l)~ and 

Except above models, there are several popular alternative 
ideas for the accelerating universe by modifying general relativ- 
ity rather than resorting to some kinds of exotic fluid. It is a nat- 
ural consequence to modify the theory of gravity, after the cur- 
rent observations mightily suggesting the failure of General Rel- 
ativity (GR) as a large cosmological scale gravity theory. Actually, 
on the right side of Einstein's field equation, the energy momen- 
tum tensor is related to the matter content, and on the left side,the 
Einstein tensor consists of pure geometrical term s, which can be 
modified to account for DM and DE phenomena dCopeland et al.1 
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2006). Accordingly, both of accelerating behaviour (DE) and dy- 
namical phenomena (DM) can be interpreted as curvature effects. 
Recently, a lot of people attempt to modify the General Relativ- 
ity in different ways and propose various modified gravity the- 
ories to address the DM or DE phenomenon or both. Amongst 
them, f(R) gravity theories, where R is curvature scalar, have 
attracted much attention, which are put forward to account for 
DE by adopting an arbitrary function of R in Einstein-Hilbert La- 
grangian, the so-called f(R) term, instead of R in traditional Gen- 
eral Relativity. Obviously, an important reason for this interest 
is that f(R) term in types of negative as well as positive pow- 
ers of R can produce the inflation at early times and now ob- 
served acceleration phase at late t imes, following the well-known 
sequence of universal evolution iNojiri & Odintsov 2003, 2007 
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Erickcek et all2006fclNayarro & Acolevenl2006tl01moi2005ll2007l 
Amendola et alj 120071 ; lli B. & Barrowl 120071 : IlI B. & Chu M-d 
2006). Usually, two possibilities referring to as the metric approach 
and the Palatini approach can derive the Einstein field equation 
from the Hilbert action. By comparison, the metric approach as- 
sumes that the equation has only one variable with respect to metric 
where the affine connection is the function of metric, while in the 
Palatini approach the metric and the connection are treated to be 
independent of each other as two variables. When f(R) is in form 
of linear function, both approaches lead to the same results. Once 
adding any non-linear term to the Hilbert action, however, the two 
methods would produce enormous differences. The reasons for the 
Palatini approach seeming appealing compared with another one 
are that, on one hand, the metric approach leads to the fourth order 
field equations, which is difficult to solve analytically, whereas the 
Palatini approach results in the two order field equations. On the 
other hand, although f(R) theories via metric approach can give 
some interesting and successful results, some of them suffer from 
certain fatal defects. For example, they cannot pass the solar system 
test, have instabilities, have incorrect Newtonian limit and cannot 
totally describe all epoches of the universe. 

Here we adopt the f(R) gravity theory within the Palatini 
approach which can avoid the above mentioned problems. In this 
framework, the form f(R) = R + a(—R) /3 is chosen so that it 
can pass the sol ar system test and has the correct Newtonian limit 
( Sotiriou l2006d), a nd can explain the late accelerating phase in th e 
universe dFav et alj2007llCapozziello et alJl2005hlSotirioull2006iJ) . 
Furthermore, it is important to go beyond the quantitative studying 
and test these theories using the observations. Recently constraint 
from data on above type of f(R) theory is inte nsively discussed by 
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2008; 



2006 



Sotirio u 2006a), the CMB shift parameter, supernovae la surveys 



data and baryon acoutic oscillations were combined to constrain 
the param eters and , in pa rticular, they give (3 ~ 10" 1 . In a previ- 
ous work. lKoivistol d2006l) used the matter power spectrum from the 
SDSS to get further restriction j3 ~ 10~ 5 , which reduced allowed 
parameter space t o a tiny around the ACDM cosmology. Recently, 
iLi B. et al] d2007t) jointed the WMAP, supernovae Legacy Survey 
(SNLS) data and Sloan Digital Sky Survey (SDSS) data to tighten 
the parameter up to [3 ~ 10~ 6 , which made this model hard to dis- 
tinguish from the standard one, where the parameter (3 is zero. So 
far, there has been no attempt to use the strong lensing observation 
to test the f(R) gravity theories in the Palatini formalism. 

Nowadays, with more and more lens surveys available to en- 
large the statistical lensing samples, gravitational lensing has de- 
veloped into a powerful tool to study a host of important subjects 



on different scales in astrophysics, from stars to galaxies and clus- 
ters, further, to the large structure of the universe. Since the grav- 
itational lensing phenomena involve the source information, two 
dimensional mass distribution of the lens and the geometry of the 
universe, it is useful to not only infer the distant source properties 
far below the resolution limit or sensitivity limit of current obser- 
vations and offer an ideal way to probe the mass distribution of the 
universe, but also constrain the p arameters of cosmological models 
dCopeland et al.l2 006; Wu 1996). So far the large lens surveys have 
built both radio (the Cosmic-Lens All Sky Survey, CLASS) and op- 
tical (SQLS DR3) lens samples. The CLASS forms a well-defined 
lens sample containing 13 lenses from 8958 radio sources with im- 

age separations of 0.3 < 8 < 3 and the i-band flux ratio q < 10 

(bright to faint) dBrowne et ai]|2002h . Meanwhile, the SQLS con- 
structs a statistical sample of 1 1 lensed quasars from 22683 optical 
quasars, of which the range of redshift is 0.6 < z < 2.2 and th e 
apparent magnitude is brighter than i — 19.1 dlnada et al]|2008l) . 
As compared with CLASS, the latter sample has larger image sep- 
arations (1 < 8 < 20 ) and smaller flux ratio limit (q c = 10~ 0,5 , 
faint to bright). Many authors use the strong gravitational lens- 
ing statistics to s t udy DE, including the equatio n of state of DE 
dChae et alJ[20o3: IChenl [20041 : lOguri et alj|2008l) a nd the test for 
the m odified gravity theories as alternatives to DE dZhu & Serend 
120081) . 

The main goal of this paper is to investigate the constraints 
of strong lensing observation on the parameters of f(R) theory in 
the Palatini approach using the SQLS statistical sample. But there 
has a big question that how the non-linear terms in f(R) gravity 
theories contribute to the light-bending and correct the Einstein de- 
flection angle. If this influence is large enough, it is inevitable that 
all the models built on GR to compute the lensing statistical prob- 
ability must be changed. But this question is starting to be solved, 
here, we still use the same way as that based on GR to calculate the 
lensing probability because of the following considerations. Firstly, 
the f(R) gravity theories are applied to explain the nature of DE, 
but as far as we know it, DE smoothly fills the space to power the 
acceleration expansion as the same wa y in everywhere of the uni- 
verse and make the universe to be flat. iDave et al, showed 
that DE does not cluster on scales less than lOOMpc. Naturally, as 
an alternative to DE, any modified gravity theory, like f(R) non- 
linear term, should not affect the gravitational potential distribution 
of dark halos as the virialized systems such as galaxies and clus- 
ters, and in particular, should have not detectable contribution to 
the light-bending if they are some successful models. Moreover, 
it is worth stressing that how the modifying the lagrangian of the 
gravitational field affects the standard theory of lensing (built on 
GR) is not well investigated. Recently, theories with f(R) oc R n 
were been investigated within the metric approach for a point-like 
lens and results were given that the R n modified gravity signatures 
possibly to b e detected through a car eful e xamination of galactic 
microlensing (Capozziello et al. 2006). And lZhangl d2007l) derived 
the complete set of the linearized field equations of the two New- 
tonian potentials <f) and ip in the metric formalism, and predicted 
that, for some f(R) gravity models, the gravitational lensing was 
virtually identical to that based on GR, under the environments of 
galaxies and clu sters. While, under the Palatini approach, the paper 
(Ruagiero 2008) is a first step to evaluate some effects of the non 
linearity of the gravity Lagrangian on lensing phenomenology. As 
expected, for a spherically symmetrical point-like model, the esti- 
mated results suggest that the effects of f(R) are confined around 
a cosmological scale, and hence, they have no effects on smaller 
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scales such as our Galaxy. For galaxies, the weak-field approxima- 
tion is valid. The standard lensing theory (based on GR) is on the 
basis of this assumption. In GR, the deflection angle for a point 
mass can be easily obtained. For the weak-field assumption, the 
field equations of GR can be linearized. Therefore, we can first di- 
vide the whole galaxy into small elements, and each mass element 
can be treated as point mass, then the deflection angle of the general 
galactic mass distribution models can be obtained by the sum (or 
integr ation) of the deflection s due to the individual mass compo- 
nents jKochanek et"alll2004h . In f(R) gravity, if the gravitational 
field is weak, we are able to perturb the General Einstein field equa- 
tions (based on f(R) gravity theories) and simplify them to linear 
equations. In this case, the superimposition principle works, so we 
can get that the field equations which describe the point-like objects 
can also be used to describe the ensemble of mass points. That is 
to say, both of the mass point and the ensemble of mass points 
have field equations in the same formalism. So we can safely ex- 
trapolate the conclusion for point-like lenses, that standard lensing 
theory is a good approximation in the case of f(R) gravity theo- 
ries, to galaxy-scale lenses. However, there are no works to prove 
it through mathematical process, and further theoretical studies are 
needed to test the f(R) gravity theories in the Palatini formalism. 
Under these situations, we consider that the standard theory of lens- 
ing is still correct in the caseof f(R) theory we used. 

According to astronomical observations, dynamical analysis 
and numerical simulations, there are three kinds of popular mass 
density profile of dark halos as a lens model in the standard lensing 
theory, namely, the singular isothermal sphere (SIS), the Navarro- 
Frenk- White (NFW) and the generalized NFW (GNFW) profile. 
The y can reasonably rep r oduce the results of str ong lensing sur- 
vey dLi & Ostrike3l20oi [20031 ; ISarbu et alj|200lh . As been well 
known, the lensing probability is very sensitive to the density pro- 
file of lenses. Moreover, the lensing is determined almost entirely 
by the fraction of the halo mass that is contained within a fiducial 
radius that is ~ %4 of the virial radius. Since in the inner regions 
the slope of SIS density is steeper than that of NFW density, for 
the multiple image separations of A9 < 5", the lensing probability 
for NFW halos is lower than the corres ponding probability for SIS 
halos by about 3 order s of magnitude dLi & Ostriker 2002,, 1200 3l ; 
iKeeton & Madaull200ll; IWvifhe et ai]|200ll) . The cooling mass di- 
viding galaxies and clusters for transition from SIS to NFW is 
M c ~ 10 13 M Q . The SQLS DR3 has 11 multi-imaged lens sys- 
tems, including ten galaxy-scale lense s and one cluster- scale lens 
with large image separation (14.62") (Inada etal. 2008). It is the 
first lens sample that contains both galaxy-scale lenses and cluster- 
scale lens. In lensing statistics, if we consider lensing halos with all 
mass scales, from galaxies to clusters, t hen the effects of substruc- 
tures dOguril 120061) and baryon infall dKochanek & White! I200l1 ; 
iKeetonlfeOQll) cannot be neglected, in particular to galaxy groups, 
for which neither SIS nor NFW can be applied. For early-type 
galaxies, which dominate the galactic strong lensing, SIS is a good 
model of lenses, in particular, it results f rom the original NFW-like 
halos through the bary on infall effects dRusin & Kochanekll2005l ; 
Koop mans et al.| [2006). Moreover, the non-spherical l ens profiles 
do not significantly affect the lensing cro ss section dOguri et au 
l2005l : lHuterer et al.ll2005l ; lKochaneki ri996). Therefore, in this pa- 
per, we use the SIS model to compute the strong lensing probability 
for galaxies. To avoid the complexity from large mass scales men- 
tioned above, we consider only galaxy-scale lenses in the sample, 
and thus compute the differential lensing probability rather than the 
usual integrated lensing probability to match the SQLS DR3 sam- 
ple (containing 10 galaxy-scale lenses). A more realistic lensing 



model including SIS and NFW to fit the complete SQLS sample 
and compari ng to the re s ults o f this paper is our future work. As 
comparison, lOguri et "all J2008h adds two additional cuts to select 
a 7 lensed quasar sub-sample from the DR3 statistical lens sample 
for constraining the parameter uo of state of dark energy. 

The rest of the paper is organized as follows: In section 2 we 
outline the f(R) theories in the Palatini approach, including the 
cosmological dynamical equation and parameters with the FRW 
universe setting. In section 3 we present the differential lensing 
probability based on the f(R) gravity theories with SIS model. 
In section 4, we discuss the specific analytic function f(R) = 
R — o>Hq(— tj^-)' 3 used in this paper, give the corresponding cos- 

mological evolution behaviour and investigate the observational 
constraint on our f(R) theory arising from SQLS DR3 statistical 
lens sample. Finally, we give discussion and conclusions in section 
5. 



2 THE F{R) GRAVITY THEORIES AND 

CORRESPONDING COSMOLOGY PARAMETERS 

2.1 The generalized Einstein Equations 

The starting point of the f(R) theories in the Palatini approach is 
Einstein-Hilbert action, whose formalism is: 



(1) 



where k — 8ttG, light velocity c = 1 and Sm is the matter ac- 
tion which is a functional of metric (? M „ and the matter fields ipM ■ 
As mentioned previously, the metric and the affine connection are 
two independent quantities in the Palatini approach, so the Ricci 
scal ar is defined b y two variables, the metric and the connection 
(see lVollickll2003l for details): 

R = g^R^, (2) 
in which 

?? — r a — r a 4- f° v x — v a v x fv\ 

where R is negative. From Eq.([T), we derive the generalized Ein- 
stein equations. Varying with respect to metric g^v, we get one 
equation: 

f'(R)R^ - \g^f{R) = -kT m „ (4) 

where f'(R) = df/dR and the energy momentum tensor T M!/ is 
given as 

2 5Sm 



(5) 



V=9 5g^ ' 

Taking the trace of Eq.© gives the so-called structural equation 
and it controls the solutions of Eq.lO- 

Rf\R)-2f(R) = -KT. (6) 
Varying with respect to the connection r A M „ gives another equation 

V Q [/'(i?)^3 M 1 = 0, (7) 

where V denotes the covariant derivative with respect to the affine 
connection. From this equation it is found that the new metric h^ v , 
which can describe the affine connection as the Levi-Civita connec- 
tions, is conformal to 



V- = f'(R)g f _ lv . 



(8) 
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Finally combined above equations, we could easily get the relation 
between R^ v and which is the Ricci tensor associated with the 
metric g M „. 

n _ n 3 V„/'V„/' V M V,f 1 vv^r 

Note that the covariant derivative V we refer to is associated with 
the Levi-Civita connection of metric g^ v . 



2.2 The Background Cosmology 

In this subsection, we shall make a detailed study of the cosmo- 
logical viability of the model based on f(R) theories in Palatini 
formalism. It is well-known that WMAP's data is a strong evi- 
dence to support a flat universe, so we choose a spatially flat FRW 
(Friedmann-Robertson- Walker) metric to describe our background 
universe. The metric takes the standard form: 



ds 2 = —dt 2 + a(t) 2 Sijdx* dx 3 (i, j = 1, 2, 3), 



(10) 



where a(t) is the scale factor. As usual, we assume a perfect fluid 
energy-momentum tensor T^" = diag(p,p,p,p). By contracting 
Eq.((9} we get the generalized Friedmann equation: 



(H ' J " 2 



lf\ 2 _l K(p + 3p) If 
2 /' 







6 /' 



(11) 



In this paper we consider a matter dominated universe, so the con- 
stant equation of state is p = u!p{uj — 0) and the relation between 
the matter density and scalar factor is p = po-R _3 '" +1 ' = 0). 
Using Eq.©, we can obtain: 

a(fi) = (/ t p m0 )3(ii/'-2/)-3 J (12) 

where we have chosen ao — 1. On the other hand, from Eq.lJSJl and 
the conservation of energy T m „ ; a = 0, we have as: 

3Hp M 



R 



(13) 



Rf"(R)-f'(RY 
Using Eqs.l(6j, i ll It and d 1 3 fc we get the Friedmann equation H(R) 



in the form of R, 

1 Rf - 3/ 



H 2 (R) 



6f n _ a f"w-2f) V : 

J \ L 2 f'(Rf"-f) I 



(14) 



2 f'(Rf"-f) ■ 

Combing Eqs. J12b and d 1 4b we can know the whole expansion his- 
tory that is determined by H(a), for any specific expression of 

f(R)- 

Let us now consider the cosmological distance based on the 
given model. From the explicit expression for the Hubble parameter 
H(R) and the relation between z and a(R): 1 + z = a~ 1 (R), it 
is convenient to rewrite the proper distance, luminosity distance, 
angular diameter distance and the deceleration parameter in terms 
of R. The proper distance and luminosity distance to the object at 
redshift z are 



D P (z) 



dz 



(15) 



Bo 



Rf" - f dR _ n p (m 
Rf — 2/ H(R) [ h 



D L (z) 



= (1 + *) 



dz 
o H(z) 



2/)- 



Rf" - f 



(16) 
dR 



(Rf'-2f)i H(R) 



= D L (R). 

The angular diameter distance from an object at red-shift z\ to an 
object at red-shift Z2 is 



D (zi,z 2 ) 



dz 



(17) 



_ h R f- 2 f)-i r R " Rf "- fl dR 

* W " (Rf>-2f)%H(R) 
= D A (R U R 2 ). 
The deceleration parameter is 

when q > the universe is in early-time matter-dominated phase 
and q < the universe is in dark energy dominated late-time ac- 
celeration phase. 



3 THE LENSING PROBABILITY 
3.1 the SIS lens 

As mentioned, in this paper, we assume that the f(R) theories have 
no effect on the gravitational lensing and still model the early-type 
galactic lens as SIS as in GR, but these formulae are used in terms 
of R. The density profile is 



p(r) 



(19) 



2ttG r 2 ' 

where a is the velocity dispersion and r is the distance from the 
galaxy center. From the geometrical relations between the image 
position 8, source position /3 and lens, one gets the lensing equa- 
tion: 

D A 

P = 0- -Ma. (20) 
u s 

For a SIS lens, the deflection angle is a — 4-k(<j/c) 2 , which is 
independent of the impact parameter. The angular Einstein radius 
is defined as 



(21) 



The sources are multiply imaged when < 6e and a ring-like 
image occurs when /3 = 0. When the multiple images occur, from 
the lensing equation, there are two images on opposite sides of the 
lens at angular positions 



9± = E ± p. 

The corresponding magnifications are 

_ 0± d0± _9 E ±f3± 

M± ~ Jl ~dW ~ At ' 

the total magnification of the two images is 



p-tot = p+ + p- 



26 E 

P 



(22) 



(23) 



(24) 



The magnification of an image is defined by the ratio between the 
solid angles of the image and the source, namely, the flux density 
ratio between the image and the source. Then the flux density ratio 
q for bright-to-faint is the ratio of the corresponding absolute values 
of the magnifications 
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<? = 



— 

/'- 



+ /3 



(25) 



Every lens survey has its own spectrum range and selection func- 
tions, so there exists an explicit cut flux density ratio, less than 
which, the survey couldnot identify the lenses. For example, for 
SQLS, the ratio (faint to bright) is larger than q c = 10" ' 5 , which 
means that only sources with ft < ftmax < 9e can be detected, 
where f3 m ax corresponds to q c . The lensing cross section is 



asis(a) = 16-7T ( — ) ( 



z f a s a,D ls D l 2 



D A 



and the velocity to produce image separation A8 is 

(A8")K 



a A e = 4.392 x 1(T 4 (— )( 



X. 



(26) 



(27) 



3.2 the differential probability 

The differential probability for the source quasar at the redshift z s 
(corresponding to R s ) lensed by foreground dark halos with mul- 
tiple image separations AO and flux density ratio larger than q c is 
given by 



P(A6,R s ,>q c ) = 



dP(> A8,R s ,>q c 



dA6 



(28) 



dR 



dD' 



dR 



^'n(a, R z )a- 3 (R)]a real (a, R z 



da 
~dA6 



where D p is the proper distance from the observer to the lens. 
n(a, Rz) and n(a, R z )a~ 3 (R) are the comoving number density 
and the physical number density of galaxies at R z (corresponding 
to redshift z) with dispersion between a and a + da, respectively, 
the cross section a rea i(a, R z ,> q c ) is dependent on the flux den- 
sity ratio of multiple images (q c ) and take into account the magni- 
fication bias (see below for details). 

In the statistics for gravitational lensing, there are two inde- 
pendent ways to get the mass function of virialized halos. One 
is the generalized Press-Schechter (PS) theory, and the other is 
Schechter luminosity function. Since the PS theory in the frame- 
work of f(R) theories does not exist, we adopt the latter. It is es - 
tablished that dMaoz & Rixll 19931 : iMoller et alj|2006l IChaejl2007h . 
the early-type galaxies dominant strong lensing, the contribution 
of the late-type galaxies is well neglected in particular when the 
image separations are larger than 1", the reason is that the late- 
type galaxies have larger rotation velocity while have smaller dis- 
persion velocity, which in turn, lead to smaller mean image separa- 
tions. Moreover, in our used SQLS lens sample all lens galaxies are 
early-type. Many previous studies for lensing statistics applying the 
no-evolution of the velocity funct ion get appealin g results to agree 
with the galaxy number co unts dim et al J 12003) and the redshif t 
distribution of lens galaxie s Jchae & Madl2003l:IOfek et al.ll2003l). 
Meanwhile, other stu dies dMao & Kochaneklll994l ; iMitchell et al.l 
l2005l :l Rix et al ] fl994h have investigated the effects of evolution of 
the velocity function on the lensing statistics, and they concluded 
that the simple evolution does not significantly affect lensing statis- 
tics if all galaxies are early-type. In our paper, we do not take 
into account the evolution of the velocity function, instead, we use 
the modified Schechter function daiaell2007l: Ichen & Zhad|2006l : 
lOguri et al.ll2008l: IMitchell et al.f2005; Zhu & Sereno 200i) 



(a)da = M~) 
a t 



' exp 



a * 



P da 



(29) 



with 



,a,p) 



r(a//3) a ' 
(0.008h 3 Mpc -3 , 161kms _1 , 2.32, 2.67) 




Figure 1. The cosmology parameter as a function of R. The upper one is 
the hubble parameter H(—R). The other one is the scale factor a(— R). 



which is derived bv lChoietal.ld2007h upon the latest much larger 
SDSS DR3 data. Then the comoving number density is n(a, R z ) = 
4>{a), and a is in terms of R z . 

Lensing probabilities must be corrected with the magnification 
bias, which explains the fact that intrinsically quasar sources with 
flux below the flux-limited survey can appear by virtue of lensing 
magnification. The cross section a rea i [a, R z ,> q c ) includi ng the 
magnification bias parameter B is defined as (Og uri et alj2008l) 



(III 



4>{Lmin/y) 
/J.(j>(L m in) 



(30) 



here du — D^pdpduj(uj G (0, 27r)) is the area element of the 
lens plane. The m agnification ratio /i in the best-fitted form is 
dOgurietal.f2 006) 



/i = MMtoi + (1 — 

here 

Jx = -[1 + tanh(1.76 - 1.780") 



(3D 



(32) 



By defining t — (3/9e and using Eqs.J2U. ( 1251 ), l |26l > and ( 1301 ), we 
can rewrite the cross section as, 

a r eai(a, R z , > q c ) = a S is(Rs,A8)B(R 3 , AO, q c ), (33) 
here the bias B for the multiply imaged source at scalar R s is 

tdt <t>{L m in/n) 



B{R s ,A8,q c 



,(34) 
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where tc = (1 — ? c )/(l + 9c) and <f>(L m in) is the quasar differ- 
ential luminosity functions at the luminosity limit L m in of the sur- 
vey, which is described by a typically doubl e power-law in terms 
of absolute magn itude in g-band at R z by ^Hopkins et al.ll2007l ; 
iLi G.L et alj|2007» 



^ M b) 1Q OA(l+0 h )(M 3 -M*) + 10 0.4(l+/3 ft )(Af g -M 3 *) ' 



and 



M* g {R s 



= MZ(Ro)-2.5[ki(a- 1 (R s ) 
+k 2 {a- 1 (R s )-l) 2 ]. 



1) 



(35) 



(36) 



For the low redshift (z < 2.1) quasars, the pa- 
rameters from the 2dF-SDSS Luminous Red Galaxy 
(LRG) and Quasi-stellar Object (QSO) Survey fitted in 
the SDSS g-band are (0*, l3 h ,Pi, Mg(Ro), fa, fe) 

-3.31,-1.45,-21.61 + 



" z mag 1 



(1.83 x 10- 6 {h/0.7) 3 Mpc 
5log(h/0.7), 1.39, -0.29). By virtual of the relation between 
luminosity and absolute magnitude L oc iq-° a ( m -m ) j we 
have (f>(L/[i) in terms of absolute magnitude. The SDSS quasar 
survey have flux limit in i-band and the apparent magnitude limit 
is imax = 19.1, but the parameters in above luminosity Equations 
are given in g-band, so we first convert the apparent magnitude to 
the absolute magnitude, then compute the corresponding absolute 
magnitude in g-band using the K-corrections which is given by 
jRichards et"aU2006h 



M g (Ro) = Mi(R ) + 2.5a, log 



( 4670A \ 
\ 7471 A J 



0.187, 



(37) 



where ct v — —0.5. 

From the differential lensing probability, we calculate the ex- 
pected number of lensed quasars from the SQLS sample, and com- 
pare to the observational results in order to constrain the parameters 
in f(R) models. The results will be given in the next section. 



4 NUMERICAL RESULTS AND OBSERVATIONAL 
CONSTRAINTS 

Here, what we would like to investigate is what extent observa- 
tions of strong lensing allow deviations from the GR which is cor- 
responding to f(R) = R. For our purposes, we ad opt the following 
representation for f(R) dAmarzguioui et alj2006l) . 



f(R)=R 



aH%{- 



R_ , 

HI' 



(38) 



which is a possible candidate for the late-time cosmic accelerating 
expansion found recently. Here, Hq is a constant with dimension 
Mpc" 1 , which are introduced to make the a and [3 dimensionless. 
Not all combinations of a and are agreement with a flat universe 
with matter dominated era flowed by an accelerated expansion to- 
day. We now consider the constraint on these parameters. Firstly, 
at the early time of matter dominated universe, the dark energy is 
ineffective and the universe is better described by GR. Therefore, 
with — R (R is negative ) being more and more larger, the modified 
Lagrangian should approach its GR limit, in order to avoid con- 
fliction with early-time physics such as Big Bang Nucleosynthesis 
(BBN) and CMB, and hence it is straightforward to demand that 
(3 < 1. Moreover, in the case of T — — p, we must demand that 
the left side of Eq.((6]l and the right of Eq.J14t be always positive. 
In the special case of (a, (5) — (—4.38, 0), the ACDM cosmology 
model is recovered. 




Figure 2. The background dynamical behaviour based on the f(R) = 
R — a(—R)@ model. The left panel is the hubble parameter H as a 
function of scale factor a. The right panel is the deceleration parame- 
ter varies with a. Different kinds of lines represent a series of values of 
/3 = -0.5, -0.1, 0, 0.1, 0.5 with n m0 = 0.27. 




A0(arc second) 

Figure 3. The expected number density based on f(R) = R — a(—R)@ 
for different choices of (3. The upper panel is for the cases of p = 
-0.005,0,0.005 and the lower panel for f3 = -0.1,-0.05,0,0.05,0.1 
from top to down. In all cases, O m o = 0.27. The case fi = corresponds 
to A CDM model. 
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ACDM 

a— 4.193,8-0 

</„> :.'fi 



A0(orc second) 

Figure 4. The expected number density in the present model for ACDM, 
our best-fitted standard model with the value of (a, /3) = (—4.193,0) 
and the ACDM in the case of (z a ) = 1.427 respectively. The two former 
number density lines based on the bin-method. 



4.1 the background evolution 

Before using Eqs.lll2t and J14t to determine the cosmological evo- 
lution behaviour, we must give the initial conditions: (p m o, Ho, 
Ro), it is noting that a subscript denotes evaluation at the present 
time. By putting /(-Ro) into Eq.d 14t. we find that the combination 
Ho/Ro appears as a single quantity in the equation, in turn, we 
choose Hq = 1 to compute Ro. Hence, from oo = 1 and the ob- 
tained value of Ro, the values of p m o and Q. m o = upmo / '(SHq ) 
are fixed. In other words, amongst the variables a, f3 and Q m o, two 
of them are independent. So far, we have obtained all the initial 
quantities. 

Now let us use the equations we derived to examine the evo- 
lution history of the universe at late-time. With the present matter 
component value of Q m o = 0.27, the changing of the Hubble pa- 
rameter and scale factor with the curvature are plotted in FigureQ] 
It is easy to see that for any choice of (5, the smaller the curvature 
\R\, the larger the scale factor, and the smaller the rate of change of 
a (a'(R)). Contrarily, the Hubble parameter H drops with decreas- 
ing \R\. Therefore, it is worth noting that the above mentioned con- 
straint conditions for further restricting the parameter (a, f3) space 
are satisfied as long as the present curvature Ro meet these condi- 
tions. And curvature Ro decreases with increasing f3. 

The background dynamical behaviour is well described by 
lines plotted in figure [2] Obviously, the Hubble parameter drop 
sharply with the expanding universe, at the same time, the value 
of deceleration parameter is changed from the positive to negative, 
representing that the universe evolves from deceleration to accel- 
eration. Different choice of ft shows different evolution history. At 
any given time, the larger value of /3 harmonizes the larger Hub- 
ble parameter and, the larger deceleration parameter after a certain 
time. It is emphasized that larger value of /? corresponds larger ac- 
celeration for our accelerated expanding universe at present time. 
Furthermore, with /3 increasing, the time for the transition from the 
gravity-dominated era to the DE-dominated era is gradually late. 



4.2 the SQLS sample constraints 

In this section, we consider the constraints arising from the strong 
lensing observations on our gravity model. We have calculated the 
number density of lensed quasars of the present model using the 




ACDM 

a = -3. 777,0=0. 0619 



A0(orc second) 

Figure 5. The lensed source quasar number as a function of image separa- 
tion A9. The histogram is the number distribution of the SQLS DR3 statisti- 
cal lens sample, in which there is a cluster-scale lens with image separation 
14.62". We rule out it because our statistical probability is computed based 
on the SIS modeled galaxy-scale lenses. The bin-size is 0.5", so the lines 
stand for N(A9) = ■ 0.5", representing the lens number around A8 
with width of 0.5". The solid line shows ACDM model, while the dashed 
line stands for our best-fitting results (a, (3) = (-3.777, 0.06195). The 
dotted line at AO = 1" indicates the SQLS resolution limit. 



SQLS sample. Specifically, from the differential lensing probabil- 
ity P(A8, Rs), which is multiplying likelihood with respect to two 
variables, image separation AO and R s (determined by the red- 
shift z s of the source quasars), if we know the number distribution 
N(z s ), then the expected number density is given by 



dN(A0)/dA0 = N(R s )P(R s ,AG)dR a 



(39) 



and the expect lens number with image separations between A8'{ 
and A#2 is computed by 

fA9 2 



N = 



dN(A6) 
dA6 



dA0, 



(40) 



which can be used to compare with 1 1 lenses in the sample to con- 
strain the parameters of the f(R) model. Furthermore, it is impor- 
tant to point out that, according to the integral mean value theorem, 
we can also rewrite the expected lens number with image separa- 
tion ranging from A#i" to A62" as follows: 



■ A0i) 



(41) 



dA6 

where AO is one exist value between A#i and A#2. 

The sample we used has 22683 quasars, and the redshifts 
range from 0.6 to 2.2. To get N(z 3 ), we count the source quasar 
number for each redshift bin, with bin-size Az s = 0.05. We in- 
vestigate how the strong lensing probability depends on the f(R) 
model we choose. In Figure [3] we present the lens number density 
based on SQLS lens sample with a series of choices of /3. When /3 
is small enough, for example /3 = ±0.005, the number density is 
hard to distinguish from ACDM (fl = 0). Inevitably, we can con- 
clude that the strong lensing statistics restricts the parameter j3 to 
not beyond the order of 10 -3 . For any choice of (3, the number den- 
sity sharply increases with small image separations, and is slowly 
down with large image separations. At about A6 = 5", the lines 
drop to zero. Different choices of the parameter significantly in- 
fluence the curves near the peak, and /3 increases while the peak 
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declines, but they cannot change the position of the peak, which 
is located approximately at AO = 0.6" and is beyond the image 
separation range of SQLS lens sample. 

We derive constraint on the parameter a assuming f) = 0. 
The best-fit value is a — —4.193, so we get f2 m o = 0.301. From 
A\ 2 = 1 we can further get the 95% confidence interval for a that 
is [-4.633, -3.754], Above approach to get the lens number den- 
sity can be called bin-method. As another interesting approach, we 
assume that the redshifts of the source quasars follow a Gaussian 
distribution, which is used to fit the redshift distribution of 22683 
quasars. The best fit gives the mean (z s ) — 1.427 and the disper- 
sion a = 0.519. In this approach, The predicted number density is 
computed by N(A0) = 22683 x P(A0, (z s ». Which is also plot- 
ted in Figure [4] in order to compare the cases of, standard ACDM 
and our obtained best-fitting standard model based on bin-method. 

As mentioned earlier, the typical type of f(R) gravity theories 
in form of f(R) = R — a(-R) 13 hav e been extensively s tudied 
by using all aspects of observa tion data. Amarzguioui et al] J20061) 
combined the SNIa data from iRiess et aL 1 20041) with the ba ryon 
acoustic oscillation length scale data lEisenstein et alj[2005t) and 
the CMB shift parameter JSpergel et alj|2003f) to constrain the pa- 
rameters of above gravity model, and found that the best-fit model 
with the value is (a,/3) = (—3.6,-0.09). Similar ly, by using 
the su pernovae data by the supernova Legacy Survey jAstier et al.l 
2006) together with the baryon acou stic oscillation peak i n the 
SDSS luminosity red galax y sample lEisenstein et al.l 120(5^) and 



. x y . 

the C MB shift parameter dSpergel et alT l2003. 2007). iFav et al.l 



J2007h constrained the parameters to a = —4.63 and f3 = 0.027 
in the best-fit case. Both works found that the ACDM model is 
well within 68.3% confiden ce level a nd the constrained parame- 
ter j3 is in 10 _ mag nitude. iKoivistd d2006h calculated the matter 
power spectrum and m atched it with the measurements of SDSS 
hegmark et al.l2004bh . The result is that the favored values of /3 are 
in order of 10~ 5 . Meanw hile, the mat t er pow er spectrum is also cal- 
culated with TT C MB bvlLi B etalJ |2007l). Where the co mbined 
data from WMAP JPage et alj|2007h Iffinshaw et al]|2007l) . SNLS 
and SDSS (Tegmark et al. 2004b) is used to restrict the model pa- 
rameters and it is found that f3 is confined in a even smaller order of 
10 -6 . From this point, we can see that the matter power spectrum is 
rather sensitive to the model parameter (3, and the model obtained 
from the observational constraints on the mater power spectrum is 
seeming indistinguishable from the standard one with (3 = 0. Here, 
we extend these studying by calculating the strong lensing proba- 
bility and using the SQLS DR3 lens sample to constrain the param- 
eters. In figure [5] we have plotted N(A8) for the ACDM model 
and our best-fitting model with a = —3.777 and /3 — 0.06195, as 
well as the histogram from the SQLS DR3 statistical sample. Ex- 
cept for both small and large image-separation ends, the non-linear 
gravity model predicts less number of lenses than those predicted 
by the ACDM model, which is linear in the form of Ricci scale R. 
The cosmological parameters fl m o = 0.285 and go = —0.544 are 
derived from the best-fit case, which are broadly in agreement with 
other measurements. We also present the contours for the joint dis- 
tribution of a and /3 in Figure [6] we can see that the ACDM model 
lies in the 68.3% confidence contour, which is consistent with the 
results of other works. Figure [7] shows the la and 2a confidence 
region of the single parameter a or /3 from the projection of the 
joint confidence lever estimated from A\ 2 = 1.0 and A^ 2 = 4.0. 
We have found that the model f(R) = J? — a{— J?)' 5 is compatible 
with the lensing observational data subject to the parameter con- 
straints [-4.67, -2.89] and [-0.078, 0.202] for a and f3, respectively, 
at the 95% confidence level . Obviously, the data coming from the 




Figure 6. The 68.7%, 95.7%, 99.7% joint confidence contour esti- 
mated from Ax 2 = 2.3,6.17 and 11.8 on the plane, arising 
from fitting the strong lensing sample. Our best fit value is (a, f3) = 
(-3.777, 0.06195), which is marked with a cross, meanwhile the ACDM 
model (a, 0) = (—4.38, 0) is marked with a star. 




Figure 7. The confidence contour plotted from Ax = 1.0, 4.0 to estimate 
the region for the single parameter a or /3. The projections on the vertical 
and horizontal axises are corresponding to la and 2a confident regions for 
the parameter (3 and a, respectively. The solid line and dotted line show the 
tangents to the ellipse representing 95% confidence level, parallelling to the 
horizontal axis and vertical axis respectively. 



SQLS DR3 lens sample can only confine the parameter /3 to order 
oflO -1 . 



5 DISCUSSION AND CONCLUSIONS 

To summarize, in this work we have derived new constraints on 
/(i?) = _R — a(-R) 13 in the Palatini formalism, using the statis- 
tics of strong gravitational lensing and the lens sample come from 
SQLS DR3. We use f(R) within Palatini approach because of it 
being free from the instabilities, passing the solar system and hav- 
ing the correct Newtonian limit. This typical type of f(R) theories 
in form of f(R) = R — a(-R) 13 have been recently put forward 
to explain the late-time expansion phase dominated by dark energy 
and have been well studied by using the cosmology measurements 
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including background universe evolution parameters, matter power 
spectrum and CMB. In our paper, we assume that the f(R) gravity 
theories have hardly detectable effect on the gravitational lensing 
phenomena of galaxy-scale halos, and then the standard lensing is 
a realistic approximation of lensing in Palatini f(R) gravity theo- 
ries. It is worth emphasizing that in the SQ1S lens sample, there is 
a cluster-scale lens, we use SIS model to calculate the differential 
lensing probability to fit the other 10 lens data to avoid destroying 
the well-defined statistical lens sample. 

Considering the FRW setting, we have shown the background 
evolution for the late expansion era of matter-dominated phase 
based on our f(R) theory. The expected lens number density aris- 
ing from the SQLS DR3 sample is not as sensitive as the matter 
power spectrum to the parameter /3. The number density distri- 
butions for different values of /3 are indistinguishable to 10 -3 in 
magnitude, while the matter power spectrum changes a lot even 
if /3 v aries only by a tiny amount, suc h as an order of magnitude 
10" 5 dLi B. et al.l2007l ; lKoivistcf 2006). Therefore, the matter spec- 
trum can constrain /? to a small region and make the model hardly 
deviation from the ACDM model. But the values of our best fit 
model, a — —3.777 and (3 — 0.06195, are not much more dif- 
ferent from the standard one with a — —4.38 and /3 = 0, im- 
plying that the model studied is not significantly preferred over 
th e standard model. The results are consist ent with the works 
of lAmarzguioui et al.l f2006h . lFav et alj J2007h and lCarvalho et al.l 
<2008h . in which the best-fitting models are (a, 0) = (-3.6,0.09), 
(a, p) = (-4.63, 0.027) and (a, (3) = (-4.7, 0.03), respectively, 
and we all found that ACDM lies in the 68.3% confidence level. 
Finally, we confine the allowed values of a and j3 to the ranges 
[-4.67,-2.89] and [-0.078, 0.202] at the 95% confidence level, and 
conclude that the SQLS DR3 statistical lens sample can only con- 
strain the parameter f3 to 10 _1 in magnitude. 

In our expected number density distribution (Figure [3}, for 
any value of /3 with the fixed Q. m o = 0.27, the peak is at about 
A9 — 0.6". But the lens sample we used, with the image- 
separation ranging from 1" to 6.17", provides no data in A6 < 1", 
this significantly restrict the constraint precision of the parameters. 
The CLASS radio lens sample has the lens image separations of 
0.3" < A6 < 3". Naturally, jointing the CLASS and SQLS sam- 
ples to constrain the model will give better results. Of course, a 
more larger lens sample from futu re surveys, (e.g., Square Kilo- 
meter Array. [Koopmans et alj|2004h , will be able to make the con- 
straints on the model more stringent. Furthermore, using the SIS 
model and NFW model together to constrain the f(R) theories 
models arising from the lensing observation data containing both 
galaxy lens and cluster lens is future work. 

Finally, we emphasize that our paper is the first try to test f(R) 
gravity theories using strong lensing statistics. The reason may be 
that, the possible deviation of gravitational lensing based on f(R) 
gravity theories from that built on GR is not well studied practi- 
cally. Our results hold on the valid hypotheses that this possible 
deviation is hardly detectable, and they are not contrary to other 
works. As mentioned in the introduction, if the f(R) theories are 
successful models to explain the nature of DE, they should have 
no important effects which can be detected on the galaxy-scale and 
cluster-scale dark halos. The reason is that DE is introduced to ex- 
plain the acceleration of the universe, which is, of course, in the 
cosmological scale. Namely, DE smoothly permeates the regions of 
smaller scales, such as galaxies and clusters, then it has no contribu- 
tions to the inhomogeneity of gravitational potential for such halos. 
For the dark matter as the dominator of the gravity field, the SIS 
and NFW lens models are well consistent with the observations. 



It is therefore necessary to test the f(R) theories through verify- 
ing that the phenomenology is not contradict those observational 
results that do agree with the predictions of the standard lensing 
theory. This critical testing will be presented in future work. 
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